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We study electrical excitation of nonlinear plasma waves in heterostructures with two-dimensional 
electron channels and with split gates, and the propagation of these waves using hydrodynamic 
equations for electron transport coupled with two-dimensional Poisson equation for self-consistent 
electric potential. The term related to electron collisions with impurities and phonons as well as the 
term associated with viscosity are included into the hydrodynamic equations. We demonstrate the 
formation of shock and soliton-like waves as a result of the evolution of strongly nonuniform initial 
electron density distribution. It is shown that the shock wave front and the shape of soliton-like 
pulses pronouncedly depend on the coefficient of viscosity, the thickness of the gate layer and the 
nonuniformity of the donor distribution along the channel. The electron collisions result in damping 
of the shock and soliton-like waves, while they do not markedly affect the thickness of the shock 
wave front. 

PACS numbers: 73.2I.-b, 73.40.-c, 73.43.Lp, 73.43.Cd 



I. INTRODUCTION 

Plasma waves, i.e., self-consistent spatiotemporal vari- 
ations of the electron density and electric field in two- 
dimensional electron gas (2DEG) channels P, 0, [H, 0, Q 
can be used in different semiconductor heterostructure 
devices. One of the most important advantages of 2DEG 
systems in comparison with 3DEG systems is the possi- 
bility to realize situations when the characteristic plasma 
frequency markedly exceeds the frequency of electron col- 
lisions with impurities and phonons. This is achieved 
using the selective doping when donors are spatially sep- 
arated from the 2DEG channel, so that the electron den- 
sity can be rather high, while the electron-impurity in- 
teraction is weakened. In 2DEG systems with highly 
conducting electrodes (gates) similar to field-effect tran- 
sistor structures, the electron density can be effectively 
controlled (in particular, significantly increased) by the 
applied voltage. Another attractive feature of the 2DEG 
system is that at realistic parameters the characteris- 
tic plasma frequency falls in the terahertz range. This 
opens up the prospects of creating novel terahertz de- 
vices, for example, sources and detectors of terahertz ra- 
diation, frequency multipliers, and so on 
recent experiments d M, El El El EI 115, 
tection of terahertz radiation in and terahertz emission 
from transistor-like 2DEG systems associated with the 
resonant excitation of plasma waves were realized (see 
also [ij], El, El, EH [H| ). Results of theoretical studies of 
plasma phenomena in 2DEG systems have been reported 
in numerous publications. However, nonlinear plasma 
phenomena in these systems are studied far less exten- 
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sively. In particular, the plasma phenomena associated 
with hydrodynamic nonlinearities in 2DEG transport 
were considered in Refs. [H, [H [H [H [H, H3, HI . The 
effect of nonlinearities associated with different contact 
effects were considered theoretically in Refs. [29l [30L l3~Lj|. 
Nevertheless, nonlinear properties of 2DEG systems can 
themselves be used in device applications. This and re- 
cent progress in experimental studies of terahertz plasma 
phenomena stimulate a significant interest in nonlinear 
plasma effects in different 2DEG systems. 

In this paper, we study nonlinear plasma phenom- 
ena associated with the electrical excitation of plasma 
waves in 2DEG systems using both analytical treatment 
and numerical modeling. We consider structures with 
a 2DEG channel supplied with side contacts and a sys- 
tem of highly conducting electrodes (split gates) which 
provides an opportunity to control the 2DEG by apply- 
ing voltage signals. The structures under consideration 
are schematically shown in Fig. 1. We use the hydro- 
dynamic model for electron transport along the channel 
coupled with the 2D Poisson equation for self-consistent 
electric potential in and around the channel. In con- 
trast to previous studies of nonlinear waves in 2DEG sys- 
tems [221, [H, HH, [25|, [26| , we take into consideration the 
2D nature of the potential distributions. This allows us 
to follow the change in nonlinear waves properties with 
the transition from 2DEG systems with the gates very 
close to the channel to those with rather remote gates 
when the plasmon spectrum essentially varies. The ef- 
fect of the spacing between 2DEG and highly conducting 
gates on the linear plasma properties we studied in many 
papers beginning with th e paper by Chaplik. [| (see, m 
particular, Refs. [H, [H [34. l35l| where the role of the 
gates and contacts on the spectra of linear plasma oscil- 
lation were considered). 

As one might expect and is shown below, the 2DEG 
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FIG. 1: Schematic view of heterostructures with 2DEG chan- 
nel and with split gate electrode configurations. 

systems with different plasmon spectra exhibit different 
behavior for nonlinear plasma waves. This is akin to the 
difference in properties of gravity surface waves in shal- 
low and deep water pools [a, 0, E2, |3f| . Apart from this, 
in our numerical modeling we focus on electron scatter- 
ing on impurities and phonons as well as on the 2DEG 
viscosity due to electron-electron collisions. 



II. EQUATIONS OF THE MODEL 



Here, v — v(t,x) and S = Y^(t,x) are the electron 
hydrodynamic velocity along and sheet density in the 
2DEG channel, ip = ip(t,x,z) is the electric potential, 
tp = ip(t,x,0) is the electric potential in the channel 
(z = 0), Erf = T.d(x) is the donor sheet concentration 
in the 2DEG channel (or near it), v is the frequency of 
the electron scattering on impurities and phonons, K is 
the coefficient of viscosity associated with the electron- 
electron scattering processes, e = |e| and m are the elec- 
tron charge and effective mass, respectively, ae is the di- 
electric constant, and 5{z) is the Dirac delta function 
describing the electron confinement in a relatively nar- 
row channel (its width is assumed to be much smaller 
than any lateral size and the distance between the 2DEG 
channel and the gate electrode). The axis x is directed 
along the channel, so that the coordinates of the side 
contacts are x = ±L/2 ~ ±L g /2, where L and L g are 
the lengths of the channel and gate (see Fig. 1), respec- 
tively {L g < L). The axis z is directed perpendicular 
to the 2DEG channel. The latter corresponds to z = 0, 
whereas the gate electrode corresponds to z = W g , where 
W g is the thickness of the gate layer. In Eq. (1), we disre- 
gard the term corresponding to the gradient of the pres- 
sure in 2DEG in comparison with the electric force (see, 
for instance, Ref. [25[). The boundary conditions follow 
from the assumption that the electron sheet densities at 
the side contacts are fixed and that the electric poten- 
tials at the electrodes (the side contacts to the channel 
and the gate sections) are given functions of time. For 
the electric potential at nonconducting surfaces (between 
the electrodes) linear approximation is used. It is also as- 
sumed that at z — > — oo, the electric field dtp/dz — > 0. 
The equations of the model under consideration account 
for the generally two-dimensional character of the poten- 
tial spatial distribution, as well as the electron transport 
in the 2DEG channel considering electron scattering and 
viscosity. 

In the case of the structures with relatively small thick- 
ness of the gate layer W g (W g <C A, where A is the char- 
acteristic length of the electron density and electric po- 
tential nonuniformities), the Poisson equation (3) can be 
replaced by the following: 



We consider the hydrodynamic equations (the Navier- 
Stokes equation and the continuity equation) and the 
Poisson equation in the following form: 
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We shall consider the following case. It is assumed 
that generally the bias dc voltage V g is applied to the 
side gate sections. Apart from this, a relatively high 
transient voltage V g (t) is applied to the central section, 
so that V g (t) = V g +AV g [@(-t)+Q(t)exp(-t/t )], where 
@(t) is the unity step function and to is the characteris- 
tic time of the gate section recharging. The recharging 
time to is determined by the capacitance, Co, of the gate 
section and the pertinent resistance, Rq, of the control 
circuit: to — RqCq. The bias gate voltages result in 
a nonuniform distribution of the electron density along 
the channel, which at t < exhibits a maximum in the 
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channel center. We shall mainly consider 2DEG systems 
with uniformly distributed donors in the channel (or in 
the gate layer slightly above the channel). The case of 
strongly nonuniform doping along the channel will be 
briefly studied as well (in Sec. V). 

The linearized version of Eqs. (1) - (3) governs the 
propagation of the plasma waves with the dispersion re- 
lation, which in the case V g = const is given by the fol- 
lowing formula: 



uj[uj + i(y + Kq 2 )] = 



4i 



W g [coth(\q\W g ) + 1}' 



(5) 



Here, 10 and q are the frequency and wavenumber of a 
linear plasma wave, 



conditions, and that the electron collisions with 
ties and phonons as well as the viscosity can be ne 
In this case, we arrive at the following equations 
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The system of Eqs. (11) - (13) can be presented as 
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is the characteristic velocity of the plasma wave in the 
gated 2DEG channel, and 
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is the ac electron sheet density. In the limits of long 
waves {qW g <C 1) and short waves UjW g ^> 1), Eq. (5) 
yields the following relationships [l], y| 
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The dispersion relations given by Eqs. (5), (8), and 
(9) are similar to those for gravity waves on the sur- 
face of a water channel. In particular, Eq. (8) corre- 
sponds to the so-called "shallow water" case, whereas 
Eq. (9) corresponds to the "deep water" case (see, for 
instance, @, 0, H3, HH). This is a consequence of the 
analogy of Eqs. (1) - (3) and the pertinent hydrodynamic 
equations for waves on the surface of a liquid under grav- 
ity force. The damping rate of the plasma wave (in the 
case v + Kq 2 <C Rew) is given by 
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III. 



RIEMANN SOLUTION AND SHOCK WAVE 
FORMATION 



Consider first Eqs. (1), (2), and (4) assuming that the 
gate layer is sufficiently thin, the length of the channel is 
sufficiently large, so that we may disregard the boundary 
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As pointed out above, Eqs. (14) and (15) are akin to 
the equations governing gravity waves in shallow water 
channels. They are also identical in form to the equations 
governing isentropic flows of a compressible gas with the 
adiabatic index 7 = 2. This indicates that Eqs. (14) and 
(15) have solutions in the form of the Ricmann waves 
which, under certain conditions, can transform into shock 
waves [H, H3] ■ One needs also to point out that the so- 
lution of Eqs. (14) and (15) should describe nonlinear 
waves with the front steepening with time until the mo- 
ment when the solution becomes three- valued. Analyzing 
Eqs. (14) and (15), we shall use a standard approach (see, 
for instance, |33,l3g|), assuming that v — w(S). Consid- 
ering the latter, Eqs. (14) and (15) become as follows: 
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Equation (18) results in 
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Considering the solution corresponding to the sign 
in Eq. (19), and introducing 



so 
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we obtain 

u(£) = 2s + const. 
Assuming that £ = Eq at v = 0, we find 



u(S) = 2(s-so) = 2so( Vj^- 1 



(21) 
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Using Eq. (22) with Eqs. (16) and (19), we arrive at the 
Riemann solution: 
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i = [w(E) + «(£)]t + C(£), 



(23) 



where s(£) is given by Eq. (20), and C(£) is determined 
by the initial conditions. Hence, the electron density as 
a function of x and t is given (in implicit form) by the 
following equation 



x = so ( 3 



2 i + C(E). 
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Using Eq. (24), we can find the wave breaking time % r , 
i.e., the time of formation of the discontinuity or the time 
of formation of a shock wave. This time is determined 
by the following conditions [36j: 
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As an example, let us consider the case when the initial 
distribution of the electron density is as follows: 
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AEq exp 
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where AEo is the magnitude of the electron density per- 
turbation and a is its characteristic length. This yields 



C(E) = a J In 




(29) 



According to Eq. (22), this distribution of the electron 
density corresponds to the following initial velocity dis- 
tribution: 



v\ t =o = 2s 
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FIG. 2: Transformation of spatial distribution o normal- 
ized electron density E/Eo in the Riemann wave governed 
by Eq. (23) with initial condition Eq. (28) and AE /E = 1. 
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FIG. 3: The same as in Fig. 2 but for initial condition given 
by Eq. (31) and AE /E = 0.5. 
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Figures 2 and 3 show the transformation of the Rie- 
mann waves governed by Eq. (23) with the initial condi- 
tions given by Eqs. (28) and (31), respectively, up to the 
moments when conditions (26) and (27) are satisfied (up 
to t = tf, r ). These moments correspond to the rightmost 
curves in Figs. 2 and 3. The electron density and the co- 
ordinate in Figs. 2 and 3 are normalized by Eo and so/v, 
respectively. The obtained spatiotemporal variations of 
E = Ti(t,x) and v = v(E(t, x)) were compared with the 
results of direct numerical solution of Eqs. (14) and (15), 
and their coincidence was confirmed. 

The solutions obtained above are related to the for- 
mation of shock waves with zero front thickness. The 
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inclusion of electron collisions and viscosity (that cor- 
respond to real systems) might affect the shape of the 
shock wave front. At a sufficiently sharp wave front, 
i.e., at a small thickness of the front 5, the viscosity 
term —Kq 2 v in Eq. (f) markedly exceeds the term as- 
sociated with the electron collisions vv. Indeed, compar- 
ing these two terms and setting q ~ <J _1 , one can find 
that the effect of viscosity dominates when S < -J K/v. 
Using Eqs. (1), (12), and (13) (valid in the case of small 
gate layer thickness) and following the standard proce- 
dure [361 ] . the shock wave thickness determined by the 
viscosity can be estimated as S ~ K/s . Substituting S 
from this equation into the above inequality, we obtain 
the condition when the viscosity prevails the collisions in 
the formation of the shock wave front: K -C Sq/v = K c . 
Assuming sq — 10 8 cm/s and v — 10 12 s^ 1 , we have 
K c ~ 10 cm 2 /s. Since the characteristic value of the 
coefficient of viscosity in real 2DEG systems is K ~ H/m 
(see, for instance, Refs. [1, 0|), where H is the reduced 
Planck constant, for the GaAs channel one obtains K ~ 
15 cm 2 /s. This implies that K <C K c and that the elec- 
tron collisions should not affect the shape of the shock 
wave front as confirmed by the results on numerical mod- 
eling in the following. 

In 2DEG systems where the gate layer is not too thin, 
the thickness of the wave front S can be determined by 
the deviation of the plasma wave dispersion relation from 
the linear one (i.e., by 2D effects leading to a difference 
in the group velocities of the harmonics with different 
wavenumber q) rather than by the viscosity. Comparing 
the contributions of these two mechanisms, one can find 
that the first mechanism dominates if K < sqW 3 = Kd 
with K d ~K c (W g v/ s ). 

Thus, in 2DEG systems with moderate gate layer 
thickness, the nonlinear waves should be fairly smooth, 
so that a soliton-like wave can form [2?], [28| rather than a 
shock wave. This is confirmed by the results of numerical 
simulations based on the general system of Eqs. (1) - 
(3) demonstrated in the next section. 



IV. SHOCK AND SOLITON-LIKE WAVES: 
NUMERICAL MODELING 

Since an analysis of nonlinear regimes governed by the 
general equations of our model, which accounts for both 
the electron collisions and viscosity, i.e. Eqs. (1) - (3), 
can not be realized analytically, we solved these equations 
numerically considering for definitcness the upper struc- 
ture in Fig. 1. The following dimcnsionlcss variables are 
used: time r = tv, coordinates £ = xu/sq and i] = zu/so, 
velocity u — v/sq, electron density a = £/£oj and elec- 
tric potential <f> = eip/ms^. In these variables, the system 
of Eqs. (1) - (3) is characterized by the following param- 
eters: L g = Lgu/so and w g = W 3 v/sq — {W g / L g )C g , 
To = tgu, and k = Kvjs\. The Galerkin spectral 
method [391 ] was used in numerical calculations. 

Figures 4-7 demonstrate the evolution of the spa- 
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FIG. 4: Evolution of spatial distributions of normalized elec- 
tron density S/Eo in the 2DEG system with normalized gate 
layer thickness w g = 10 -3 . The inset corresponds to the early 
stage when the shock wave is formed. This case is analogous 
to the propagation of surface waves in a "shallow water" chan- 
nel. 



tial distributions of the electron density (normalized by 



£n) in the structures with C„ 



10 and different nor- 



malized gate layer thicknesses (w g — 0.001 — 1) in re- 
sponse to a drastic change of the central gate section 
potential: Av g — 2exp(— 100t) (t = 10~ 2 ). It is as- 
sumed that k = 2 x 10~ 3 . These cases correspond 
to the spatial distributions of the electron density at 
the initial moment with a maximum at a; = 0, where 
(£/£d)|t=o =- 1-3-3. If s = 10 8 cm/s and v = 10 12 s" 1 , 
the results demonstrated in Figs. 4-7 correspond to the 
gate length L g = 10 /im (i.e., to a fairly long channel) and 



W g /L g 



10" 



T0 . The length of the central gate sec- 



tion is about 0.5 /im. As shown, in systems with small 
values of the parameter w g (Figs. 4 and 5), the initial 
perturbation of the electron density transforms into two 
pulses propagating from the channel center to the side 
contacts. After a short time, the front of these pulses 
becomes sharp, so the pulses turn into shock waves. The 
results obtained for the case of small w g coincide with 
the results of the calculations by Dmitriev at al. [23[ and 
Rudin at al. [24|,[25[. However, in 2DEG systems with 
larger values of the parameter w g , the propagating waves 
have the form of relatively smooth pulses with a soliton- 
like shape (Fig. 6). At sufficiently large w g , the transfor- 
mation of the initial perturbation of the electron density 
(Fig. 7) is even qualitatively different from that in the 
case of effectively gated 2DEG systems. This difference 
in the nonlinear waves propagated in 2DEG systems with 
different parameter w g can be attributed to a significant 
difference in the dispersion relations of the plasma waves 
in the cases of small and large spacings W g (small and 
large values of parameter w g ) . 

As seen from Figs. 4-7, the amplitude of the propa- 
gating pulses markedly decreases with time. This effect 
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FIG. 5: The same as in Fig. 4 but for w g — 10 
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FIG. 6: The same as in Fig. 4 but for w B = 1CT 1 . 
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FIG. 8: Snapshots of spatial distributions of normalized elec- 
tron density E/Eo with oscillatory structure behind the shock 
wave front in the system with normalized gate layer thickness 
w g = 10~ 2 and small viscosity (k = 10~ 4 ). 



is obviously associated with electron collisions. Compar- 
ing Figs. 4, 5, and 6, one can see that the thickness of the 
wave front markedly increases with increasing w g charac- 
terizing the role of the plasma wave dispersion associated 
with the 2D nature of the distributions of the electrical 
potential in the systems under consideration. 

Figure 8 shows the evolution of the electron density dis- 
tributions in a 2DEG system with parameter k substan- 
tially smaller than in Figs. 4-7. One can see that in this 
case, the electron density distributions behind the shock 
wave front exhibit pronounced oscillations. To avoid nu- 
merical calculation artifacts, the pertinent modeling was 
conducted with rather different numbers of modes in the 
Fourier transformation (varied from 2 x 10 3 to 2 x 10 4 ). 
The independence of the result of the number of the 
modes used in the calculations was confirmed. 

The oscillatory character of the electron density distri- 
butions can be attributed to the instability of the shock 
waves with relatively smooth (but nonuniform) distri- 
butions behind the shock wave front. As shown previ- 
ously [40(, the nonuniformity of 2DEG in the channel 
can markedly affect the conditions of plasma instabil- 
ity. Considering the nonuniformity of the electron density 
distribution in a propagating shock wave and assuming 
that the characteristic length of this perturbation q^ 1 is 
small in comparison with the scale of the nonuniformity 
of the electron density distribution behind the front, one 
can get the following equation for the damping rate of the 
electron density perturbations behind the front (compare 
with Eq. (10)): 



FIG. 7: The same as in Fig. 4 but for w g = 1 (thick gate 
layer - "deep water" case). 
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FIG. 9: Time dependences of the front velocity for the same 
parameters as in Figs. 4-6. 
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be satisfied for the perturbations with the spatial period 
increasing with time. Just such a behavior of the oscil- 
latory structure behind the shock wave front is seen in 
Fig. 8. 

Figure 9 shows how the velocity of the wave front (nor- 
malized by so) changes with time (normalized by 
We define the velocity of the wave front as the velocity of 
movement of the point where £ = (Yi max + So)/2. Here, 
^max is the peak value of the electron density. One can 
see that in 2DEG systems with relatively thin gate layers 
in which the formation of the shock waves occurs (corre- 
sponding to Figs. 4 and 5), the velocity of the front move- 
ment exceeds the characteristic plasma wave velocity Sq 
but this front movement decelerates with elapsed time. 
This can be attributed to the fact that, according to the 
general properties of shock waves (35[ (see also Ref. 25] ) , 
their velocity is determined by so and s — sov^V^o ( see 
Eq. (20)) with the latter value decreasing throughout the 
wave propagation (as seen in Figs. 4 and 5). However, 
the velocity of propagation of more smooth pulses (see 
curve marked by squares in Fig. 9) in the case of larger 
gate thickness is close to sq. 

The results shown in Figs. 4-8 correspond to rela- 
tively fast switch-off of the voltage applied to the central 
gate section (r = 1CP 2 ). The obtained results do not 
change when To increases up to To = 10 . However, at 
relatively slow voltage switch-off, i.e. at moderate val- 
ues of parameter t , the pattern of relaxation of electron 
density becomes markedly different. Figure 10 shows the 
relaxation of electron density in the 2DEG system with 
the same parameters as in Fig. 4 at t = 1. In this 
case, the transformation of the electron density distri- 
bution closely resembles the density peak smearing con- 
sidered previously in the framework of the conductive 
model (SEl. 



FIG. 10: Evolution of spatial distribution of normalized elec- 
tron density E/Eo in the 2DEG system with the same pa- 
rameters as in Fig. 4 but at relatively slow voltage switch-off 
(to = 1). 



In the case Imw > 0, i.e., when 



V. WAVE PROPAGATION IN NONUNIFORM 
CHANNEL: EFFECT OF "TSUNAMI" 

Consider now a 2DEG system with markedly nonuni- 
form distribution of donors = T^d(x) along the chan- 
nel, so that the quantity So given by Eq. (7) also depends 
on the coordinate x. We assume that 
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the nonuniform but relatively smooth electron density 
distribution behind the shock wave front can become un- 
stable resulting in the formation of oscillatory distribu- 
tions. This instability condition shows that the electron 
density distribution can be unstable when the gradient of 
the electron density is sufficiently large. Since the mag- 
nitude of the electron density peak decreases with wave 
propagation (due to electron collisions), and the total 
number of electrons is constant, the electron density gra- 
dient decreases as well. As a result, inequality (34) can 



if \x\ < L , 



PS- 



2\x\-L a 
2a d 



if 



> L , 



where T, d is the sheet density of donors in the channel 
center, Yi e d is the donor density near its edges, and L c i 
and ad are the characteristic lengths of the spatial distri- 
bution of donors (Ld < L). Figure 10 shows the evolution 
of the spatial distribution of the sheet charge density in 
the 2DEG channel normalized by (So), where the symbol 
(...) means the averaging over the channel length. One 
can see that two soliton-like waves propagating in oppo- 
site directions toward the side contacts arise. However, 
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s 
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FIG. 11: Evolution of spatial distribution of normalized sheet 
charge (E — Eo)/(£o) in the 2DEG system with nonuniform 
donor distribution. 

upon reaching the regions with relatively low donor den- 
sity (r ~ 0.25), these waves transform into shock waves 
with rather sharp fronts and complex structures behind 
the front. In particular, the fronts are followed by deple- 
tion regions and then by the secondary shock waves. The 
behavior pattern of the wave originated from the initial 
pulse resembles tsunami in the ocean. This is due to an 
analogy between the equations governing the 2DEG in 
the systems under consideration and the equations de- 
scribing the gravity surface waves in shallow and deep 
water pools. 



result in damping of the shock and soliton-like waves. 
However, they do not markedly affect the thickness of 
the shock wave front. Due to the fairly sharp wavefront 
of nonlinear waves in the 2DEG systems under consid- 
eration or due to their strongly oscillatory structure (as 
shown in Fig. 8), the transient charges induced by the 
propagating charges of these waves in the side contacts 
and, hence, in the external circuit might exhibit steep 
surges. Fast transient currents in such a circuit includ- 
ing an antenna might provide generation of electromag- 
netic radiation with relatively high frequencies. The het- 
erostructures with a nonuniform doping of 2DEG chan- 
nel appears to be rather promising. In the structures 
with the gate separated into several short sections (see 
the lower structure in Fig. 1) the propagating charges 
of shock waves might result in even steeper surges of the 
charges induced in these section [42j ■ Since these sections 
can serve as an antenna components, fast variations of 
the induced charges and, hence, the dipole momentum, 
can lead to emission of electromagnetic radiation. There 
are several characteristic frequencies related to the sit- 
uation: /i ~ so/L gi fi r~j So/lg, where l g is the lat- 
eral spacing between the neighboring gate sections, and 
f3 ~ sq/S. The frequencies ft, $2, and /3 can be much 
larger that the inverse recharhing time t^ 1 easily fall in 
the terahertz range at realistic values of the character- 
istics lengths due to fairly large values of the character- 
istic velocity of plasma wave so (about 10 8 cm/s). The 
mechanism of the generation of relatively high frequency 
electrical signals, which convert into electromagnetic ra- 
diation is akin to the frequency multiplication due to the 
plasma nonlinearity. Thus, the excitation of shock and 
soliton-like waves in the heterostructures under consid- 
eration by external electrical signals might be used to 
generate terahertz radiation. 



VI. CONCLUSIONS 

We studied electrical excitation of nonlinear plasma 
waves in heterostructures with 2DEG channels and split 
gates. The propagation of these waves was considered 
both analytically and using numerical modeling. The 
hydrodynamic electron transport model accounted for 
electron collisions with impurities and phonons as well 
as 2DEG viscosity. The hydrodynamic equations were 
supplemented by the 2D Poisson equation for the self- 
consistent electric potential. We demonstrated the for- 
mation of shock and soliton-like waves as a result of the 
evolution of strongly nonuniform initial electron density 
distribution. It was found that the shock wave front and 
the shape of soliton-like pulses pronouncedly depend on 
the coefficient of 2DEG viscosity, the thickness of the gate 
layer, and the nonuniformity of the donor distribution 
along the channel. In the case of the 2DEG channel with 
strongly nonuniform doping, the effect of transformation 
of relatively smooth waves into shock waves, resembling 
a tsunami effect, was observed. The electron collisions 
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